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Radiation-induced  sound  experiments  were  conducted  with  one  to 
four  10.6p  CO^  laser  pulses.  These  experiments  were  conducted  to 
demonstrate  the  ability  to  control  the  acoustic  energy  spectral  density. 
The  acoustic  signatures  were  Fourier  analyzed  to  obtain  the  energy 
.spectral  density.  The  results  show  that  the  ratio  of  the  width  of  a peak  in 
the  spectrum  to  the  fundamental  frequency  scales  approximately  as  the 
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ABSTRACT 


[ Radiation  induced  sound  experiments  were  conducted  with  one  to 

! four  10.6ii  CO2  laser  pulses.  These  experiments  were  conducted  to 

. 

I demonstrate  the  ability  to  control  the  acoustic  energy  spectral  density. 

I The  acousic  signatures  were  Fourier  analyzed  to  obtain  the  energy 

* spectral  density.  The  results  show  that  the  ratio  of  the  width  of  a peak 

■ in  the  spectrum  to  the  fundamental  frequency  scales  approximately  as  the 

I inverse  of  the  number  of  pulses.  This  scaling  does  not  hold  if  the  pulses 

; are  not  evenly  spaced  or  not  uniform  in  amplitude. 

f 

I Theoretical  modeling  formulated  scaling  laws  relating  laser 

parameters  to  acoustic  parameters  of  the  generated  signatures. 

i: 

' ' Theoretical  predictions  were  also  made  of  the  change  in  amplitude  and 


shape  of  the  broadband  signals  due  to  frequency-dependent  absorption  as 
a function  of  propagation  distance. 
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I.  INTRODUCTION 


The  physics  of  radiation  induced  sound  has  previously  been  investi- 

1 -9 

gated  both  theoretically  and  experimentally  “ , but  past  experiments  have 

15  8 

been  concerned  with  single  pulses  only  ’ , In  this  report  we  discuss 

theoretical  and  experimental  results  with  one  to  four  laser  induced  acoustic 
pulses  and  the  resultant  pulse  train. 

A major  criticism  of  prior  radiation  induced  sound  experiments  was 
that  the  source  bandwidth  was  too  great  to  be  useful.  Therefore,  this  ex- 
perimental program  was  undertaken  to  demonstrate  that  narrower-band 
sound  could  be  generated  with  multiple  laser  pulses.  The  most  flexible 
and  cost  effective  approach  was  to  simulate  only  the  first  four  pulses  of  a 
long  pulse  train  with  a laser  system  consisting  of  four  independently 
triggered  lasers.  This  laser  system  was  assembled  from  four  CO2  laser 
kits  and  one  power  supply.  The  experiments  are  discussed  in  Section  III. 

The  experimental  results  demonstrate  that  the  bandwidth  can  be  controlled 
by  the  number  of  laser  induced  acoustic  pulses. 

The  theoretical  effort,  reported  on  in  Section  II,  was  an  extension 
of  previous  work.  The  earlier  analysis,  which  was  done  for  PMS-405  of 
the  Naval  Sea  Systems  Command,  established  the  basic  scaling  laws  of 
3D  blast  wave  induced  sound.  ^ This  work  was  recently  extended  and  modi- 
fied to  predict  theoretically  the  far-field  acoustic  amplitude  as  a function  of 
9 

time.  The  far-field  signature  is  predicted  from  linear  acoustics,  and  the  pre- 
dictions were  made  to  estimate  the  change  in  the  signature  due  to  the 
frequency  dependent  absorption  coefficient  of  water.  No  dispersion  of 
sound  waves  is  used  or  appears  to  have  been  observed,  but  a broadband 

4 

signature  does  spread  out  in  time.  This  spreading  is  due  to  the  approxi- 
mately quadratic  dependence  of  the  absorption  coefficient  on  frequency: 


for  a fixed  propagation  distance  the  high-frequency  components  making 
up  the  broad  band  signature  are  attenuated  more  than  the  low-frequency 
components.  It  is  well  known  from  Fourier  analysis  th.^’  a waveform 
may  spread  in  time  as  its  spectrum  is  reduced  in  width. 


Conclusions  and  suggestions  for  future  work  are  discussed  in 
Section  IV. 


II.  THEORY 


Past  theoretical  work  on  radiation-induced  sound  treated  the 

phenomenology  in  two  limits.  First,  the  steady-state  coupling  problem 

was  analyzed.  ^ This  work  yielded  the  conclusion  that  an  induced  surface 

pressure,  p , would  be  proportional  to  the  laser-beam  irradiance,  I : 
s o 


p = C I 
s o 


where  C,  the  ratio  of  the  pressure  to  the  irradiance,  is  commonly  called 

the  coupling  coefficient.  The  work  done  by  Pirri  yielded  C 25  dyne/W  for 
5 7 2 ^ 

10  < I^  < 10  W/cm^.  In  the  second  limit  of  the  phenomenology,  it  was 

assumed  that  a steady-state  situation  would  not  occur,  and  that  blast-wave 

theory  would  be  most  appropriate  to  use.  ’ The  application  of  blast-wave 

theory  to  analyze  the  radiation-induced- sound  phenomenology  is  a method  to 

uncouple  the  optical  radiation  from  the  fluid  mechanics  of  the  water  and  water 

vapor.  In  the  model,  a short-pulse  laser  beam-pulse  laser  beam  irradiates 

water  with  a fluence  much  greater  than  the  threshold  fluence.  The  threshold 

fluence  is  defined  as  the  energy  per  unit  area  sufficient  to  "boil''  an 

infinitesimal  layer  of  water  and  is  a function  of  the  optical  depth  and  the 

properties  of  water.  This  calculated  energy  must  be  delivered  in  a time 

short  enough  so  that  heat  conduction  is  insignificant.  The  calculated  values 

of  E^  are  reproduced  from  Ref.  9 in  Table  I for  various  laser  wavelengths. 

We  were  experimentally  restricted  to  10.  6 p radiation,  for  which  the  threshold 
2 

fluence  is  3 J/cm  . When  the  laser  fluence  is  much  greater  than  the  thresh- 
old fluence,  a small  amount  of  water  is  raised  to  a temperature  and  pressure 

high  in  comparison  with  the  normal  boiling  conditions  of  100°C  and  10^ 

2 

dynes/cm  . The  water  vapor  drives  a shock  wave  away  from  the  irradiated 
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TABLE  I 


THRESHOLD  FLUENCE  FOR  VARIOUS  LASER 

LINES 

.aser 

Lines 

Laser 

Wavelen^h 

um 

Threshold 

Fluence 

J/cm'^ 

CO^ 

P (20) 

10.6 

3 

CO 

8-7  Band  P (14) 

5.  25 

11 

9-8  Band  P (10) 

5.  27 

11 

DF 

2-  1 Band  P (6) 

3.  73 

14 

2-1  Band  P (8) 

3.  80 

14 

2-1  Band  P (9) 

3.  84 

14 

HF 

2-1  Band  P (7) 

2.87 

0.  2 

2-1  Band  P (8) 

2.  91 

0.  2 
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area.  We  assume  that  for  the  duration  of  the  blast  wave,  i.  e.  , for  tens 


of  microseconds,  the  behavior  is  dominated  by  the  blast-wave  solution 
in  air  and  that  the  water  surface  is  unperturbed.  Sensible  displacements 


of  the  water  occur  tens  of  milliseconds  after  the  radiation  pulse. 


Blast-wave  theory  assumes  an  instaneous  release  of  energy  with 
n omass  addition.  In  our  case  there  can  be  significant  mass  addition,  but 
it  will  be  ignored  to  keep  the  problem  tractable.  Recently,  blast-wave 
theory  has  been  extended  to  include  a noninstantaneous  release  of  energy. 
We  use  this  formulation  while  the  laser  pulse  is  on,  and  then  match  the 
solutions  of  classical  blast  wave  theory  to  the  variable  energy  theory  at 


the  pulse  "off"  time. 


We  consider  a constant-input -power  case;  the  power  can  be 
deposited  in  a plane  for  ID  expansions,  in  a line  for  2D  expansions,  and 
at  a point  for  3D  expansions.  The  shock  radius  and  pressure  are  used  to 


model  the  radiation-induced  sound  source.  Variable-blast-wave  theory 

10 

yields  a solution  for  the  shock  radius  R , which  is: 


R = Kt 
s 


3/(a  + 3) 


where 


^o  \ a + 3 / o 


l/(a  + 3) 


J is  a constant  with  values  1,  .57  and  .4  , > nd 
o 


a = 0,  1,  2 for  ID,  2D,  and  3D  expansions,  respectively. 


a is  the  sound  speed  in  air  , 
o 


-5- 


is  the  ambient  pressure  , 


I is  the  input  power  per  unit  area  for  a = 0 > 


I is  the  input  power  per  unit  length/2  tt  for  cc  = 1,  and 


I is  the  input  power  / 4tt  for  a = 2. 


t is  time  measured  from  the  start  of  the  energy  release  (i.  e. 
from  the  time  at  which  the  accumulated  laser-applied  energy  per  unit  area 
equals  the  threshold  fluence)  and  limited  its  value,  t,  at  the  end  of  the  laser 
pulse.  The  shock  pressure,  p,  is  given  by: 


p = ^ 


(2.  3) 


where  Eq.  (2.  2)is  used  for  R and  y is  the  ratio  of  specific  heats.  The  early 

s 


time  solution  will  be  one  dimensional.  In  this  case  we  have: 


(2.4) 


The  surface  pressure  for  the  ID  case,  is  approximately  one  half  the 

shock  pressure  because  the  pressure  behind  the  shock  front  falls  off  steeply  to 


a value  approximately  one  half  the  value  given  by  Eq.  (2.4).  Thus  we  have 

- ,2/3 

1 i ^ ^ \ 

p = — — ^ p 1 — 1 0 ^ts  T . (2.5) 

ID  Y + 1 ° V Pq  / 


By  numerical  computation  with  Eq.  (2.  2)  and  an  Input  power  per  unit 


8 2 

area  of  10  W/cm  the  shock  radius  will  have  grown  to  1 cm  in  Ipisec. 


Because  the  spot  sizes  which  were  experimentally  used  were  approximately 
1x2  cm,  ID  variable-energy  blast-wave  theory  is  appropriate  to  use  during 
at  least  the  first  part  of  the  laser  pulse  duration.  However,  after  the  laser 


1 


pulse  if  off  the  expansion  is  three  dimensional,  because  the  shock  radius  has 
grown  to  be  larger  than  the  spot  dimension.  In  what  follows,  the  approxima- 
tion has  been  made  that  one-dimensional  blast-wave  theory  is  appropriate 
during  the  whole  duration  of  the  laser  pulse  and  that  3-dimensional  blast- 
wave  behavior  occurs  immediately  thereafter. 

Figure  2.  1 is  a plot  of  the  surface  pressure  versus  I,  the  available 
power  per  unit  area.  These  values  will  be  used  in  the  computation  of  a 

far-field  acoustic  signature.  Experimental  measurements  have  demonstrated 

9 2 

clean-air  breakdown  levels  of  greater  than  10  W / cm  for  10.  6 p radiation. 

7 

Dirty-air  breakdown  levels  with  a 10.  6 (i  laser  wavelength  range  from  10  to 
9 2 

10  W/cm  for  particles  ranging  from  30  to  1 pi  diameter  for  a 1 i-isec 

11  8 2 

pulse,  and  an  air/water  breakdown  threshold  of  10  W/cm  with  10.  6 (i 

2 

laser  wavelength  has  been  measured. 

Blast-wave  theory  requires  specification  of  the  power  or  energy 
deposited  in  the  gas.  An  approximate  method  to  relate  the  laser  energy  to 
the  blast-wave  energy  per  unit  area  is  to  subtract  from  the  pulse  duration 
a time  sufficient  for  the  deposition  of  the  threshold  fluence.  After  the 
laser  pulse  is  off,  the  blast  wave  will  decay  according  to  the  3D  blast-wave 
theory  discussed  above.  The  shock  radius,  R^,  is  given  by 


(2.6) 


where  is  taken  to  be  the  shock  radius  at  the  end  of  the  laser  pulse 
which  has  a pulse  time  t;  Rj  can  be  evaluated  from  Eq.  (2.2).  The 
surface  pressure,  , scales  according  to 


Pg  = P (T)  (-  ) 


(2.  7) 


where  p (t)  is  the  pressure  at  time  t.  We  match  the  ID  and  3D  pressures 
at  this  time,  and  find  that 


P|D 


(2.8) 


The  above  information  will  now  be  used  to  compute  an  acoustic  signature. 
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The  well  known  formula  for  a broadband  acoustic  signature, 

p (r,  t + measured  at  a point  normal  to  the  source  is: 
a 


p (r,  t + -) 
a 


— - / 
r a dt  J 


P dA  , 
s 


(2.  9) 


where  a is  the  sound  speed  in  water  and  the  notation  indicates  an  area 
integration  over  the  water  surface. 

Equation  (2.  9)  will  be  solved  by  inserting  Eqs.  (2.  5)  and  (2.  8), 
the  assumption  that  the  surface  pressure  is  constant  over  an  area  of 
the  water  surface  and  zero  elsewhere  will  be  used  to  simplify  the  integrals. 
Thus  we  have: 


p (r,  t + “ ) - A 

^ a 2TT  ra  o 


V+  1 


6(t.-), 


0 < t < T 


(2.  10) 


p (r,  t + -)  = 

ct 


- (-)  — 
[5/  2Ttr 


~ ttR, 
a 1 


(2.  11) 


where  the  approximation  has  been  made  that  while  the  laser  pulse  is  on, 

the  area  A is  the  laser  beam  area.  After  the  laser  pulse  is  off,  the  3D 
o 

expansion  of  the  shock  is  used  over  the  water.  We  note  that  the  problem 

7 

for  a solely  3D  expansion  has  already  been  treated.  The  approximation 

which  "matches”  the  blast-wave  behavior  during  the  time  when  the  laser 

pulse  is  on  to  the  subsequent  behavior  (when  the  laser  pulse  is  off),  requires 

2 

the  assumption  that  TT  Rj  in  Eq.  (2.  1 1)  is  equal  to  A^  in  Eq.  (2.10).  A 
more  exact  treatment  is  beyond  the  scope  of  this  work. 


At  this  point  it  is  appropriate  to  consider  the  validity  of  Eq.  (2.  9), 

which  assumes  that,  in  the  water,  linear  acoustics  is  applicable.  From  the 

10  2 

bulk  compressibility  modulus  of  water  (2  x 10  dyne/cm  ) --  an  overpressure 
9 2 

of  10  dynes/cm  would  be  required  to  generate  a density  change  of  5%.  Thus 

9 

it  would  appear  Eq.  (2.  11)  is  valid  as  long  as  the  pressure  is  below  10 
2 

dynes/cm  . We  shall  see  that  this  criterion  is  met  for  our  experiments. 

We  note  that  for  a fixed  pulse  time  T,  the  time  for  the  signature 

to  decay  to  any  given  fraction  of  its  greatest  amplitude  is  independent  of 

energy.  Therefore,  only  the  magnitude  of  the  energy  spectral  density  will 

change,  and  not  the  peak  frequency  or  other  aspect  of  spectral  shape.  This 

conclusion  is  valid  only  within  the  context  of  our  idealized  model.  In  fact, 

a blast  wave  with  counterpressure  will  yield  an  energy  spectral  density 

changing  in  magnitude  and  shape  as  a function  of  energy.  However,  the 

solution  of  the  blast  wave  problem  with  counterpressure  must  be  treated 

numerically,  and  we  will  restrict  our  attention  to  the  analytic  solutions 

of  blast  waves  without  counterpressure  even  though  certain  anomalies  in 

the  resulting  solution  could  be  developed  there  from. 

A graphical  representation  of  the  far-field  term  of  Eq.  (2.  11)  is 

useful  for  making  predictions  of  the  far-field  acoustic  amplitude.  We 

2 

assume  a 1-^isec  laser  pulse  and  1-cm  laser  beam  area  and  evaluate  the 
far-field  amplitude  at  t = Ij-isec,  t = 5|isec,  and  t = lO^isec  for  r = 1 m. 

The  results  are  shown  in  Fig.  2.  2 where  the  rare  faction  amplitude  is 

plotted  versus  irradiance.  As  noted  above,  the  threshold  fluence  is 

Z 6 Z 

3 J/cm  , and  the  plot  cannot  be  used  below  3 x 10  W/cm  if  t is  Ipsec. 

The  hydrophone  may  be  considered  to  be  in  the  far  field  if  its  distance 

from  the  dipole  source  is  greater  than  one  wavelength  of  sound  for  each 

frequency  of  interest.  This  criterion  is  met  for  frequencies  greater  than 

1500  Hz.  Therefore,  for  measurement  depths  of  1 m there  is  no  doubt 

that  the  data  is  far  field  for  our  broadband  signatures. 

The  predictions  of  the  propagation  characteristics  of  the  generated 

acoustic  pulses  is  a nontrivial  problem.  An  approximate  method  will  be 

used  to  predict  the  signature  for  a propagation  path  for  which  absorption 

-10- 
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is  dominant.  The  equations  describing  the  absorption  coefficient  as  a 

12 

function  of  frequency  were  taken  from  Urick.  Two  graphs  are  repro- 
duced from  Urick  in  Fig.  2.  3,  depicting  the  absorption  coefficient  for 
40°F  water.  There  is  no  dispersion  in  the  speed  of  sound  in  pure  water. 
Because  the  dispersion  relationship  is  not  known  to  us  for  sea  water,  the 
nondispersive  nature  of  pure  water  will  be  assumed  in  our  analysis. 

The  computational  procedure  is  as  follows;  first,  the  acoustic  wave 
form  is  digitized  at  equally  spaced  intervals,  and  the  digitized  waveform  is 
plotted.  Second,  the  Fourier  transform  is  calculated  by  fast  Fourier  trans- 
form techniques.  Third,  we  compute  an  energy  spectrum  and  plot  density, 
which  is  the  squared  modulus  of  the  Fourier  transform.  Fourth,  the  Fourier 
transform  is  attenuated  at  each  frequency  by  the  amount  that  a monochromatic 
wave  would  be  attenuated  for  the  specified  propagation  distance.  It  is  this  step 
which  allowed  an  approximate  solution.  It  is  not  obvious  that  the  attenuation 
of  a broadband  signal  should  be  calculated  from  the  superposition  of  mono- 
chromatic attenuated  values.  Also,  the  assumption  that  there  is  no  dispersion 

in  the  presence  of  frequency -dependent  attenuation  is  not  realistic  for  sea  water 

13 

but  true  for  pure  water.  However,  the  lack  of  a dispersion  information 
necessitates  the  use  of  this  method.  Fifth,  an  inverse  Fourier  transform  is 
performed  to  yield  an  attenuated  signature.  The  signatures  can  be  computed 
for  arbitrary  propagation  paths.  Spherical  spreading  losses  are  not  included 
in  the  normalization  constants. 

An  example  of  the  result  of  this  procedure  is  shown  in  Fig.  2.4 
for  Eq.  (2.  10)  and  (2.  11).  The  delta  function  is  approximated  by  a 
single  point  and  requiring  that  the  integral  of  the  total  signature  be 
zero  yields  the  correct  amplitudes.  This  requirement  is  obvious  from 
Eq.  (2.  9)  since  the  time  integral  of  the  right  side  of  Eq.  (2.  9)  is  zero.  The 
propagation  path  of  1 m was  chosen  to  match  the  experimental  distance, 
and  the  example  was  done  for  T = Ipsec.  The  energy  spectral  density  is 
shown  in  Fig.  2.  5.  Another  example  was  computed  for  a pulse  duration  of 


J 


Fig.  2.3  Absorption  Coefficients  versus  Frequency  from  Urick 
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3 ^sec  and  is  shown  in  Figs.  26  and  27.  These  two  examples  were  done 
because  they  bound  the  laser  time  regime  of  the  experiments.  Signatures 
at  SOO,  1000,  1500  m are  also  given.  Spherical  spreading  is  not  included 
in  the  amplitude -normalization  factor. 

So  far  we  have  restricted  the  discussion  to  single-pulse  signatures. 

Now  we  use  a theorem  from  Fourier  analysis  to  extend  the  discussion  t<j 
multiple  equally  spaced  identical  pulses.  The  theorem  is  that  the  energy 
spectral  density  of  N pulses,  (f)  is  given  by 

sin^  N TT  f/f 

f (f)  = r 2_  * (f)  ^ (2.13) 

sin  (nf/f  ) 
o 

; 

where  f^  is  the  repetition  frequency  of  the  pulses  and  5^  (f)  is  the  energy 

spectral  density  of  one  pulse.  Thus  the  theoretical  work  of  this  section  . 

2 2 ' 

is  easily  extended  to  multiple  pulses  by  the  function  sin  Nrrf/f^/sin  nf/f  . 

An  illustration  of  this  for  one,  two,  three,  and  four  cycles  of  a sine  wave 
separated  by  a space  equal  to  four  periods  is  shown  in  Fig.  2.  8.  The 
enhanced  amplitude  and  bandwidth  narrowing  of  the  fundamental  and 
harmonics  is  apparent. 
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ni.  EXPERIMENT 


The  experimental  program  was  undertaken  to  demonstrate  the  effect 


of  the  number  of  acoustic  pulses  on  the  energy  spectral  density.  As  dis- 
cussed in  Section  II,  if  the  pulses  are  identical  and  evenly  spaced,  the 
resultant  energy  spectral  density  for  N pulses  ? (f)  will  be,  (Eq.  (2.  13); 


sin^  N tt  f/f 


Sin  rr  i/i 


(f) 


(3.  1) 


Before  the  start  of  this  program,  it  was  not  clear  whether  identically  shaped 
acoustic  waveforms  could  be  generated  with  a sequence  of  identical  laser 
pulses  irradiating  the  same  spot.  Reasons  for  not  expecting  identical 
pulses  are  discussed  below.  As  will  be  shown  in  this  section,  the  wave- 
forms that  were  generated  were  generally  very  similar  in  shape,  but  not 
necessarily  in  amplitude.  At  this  point  it  is  not  clear  whether  the  change 
in  amplitude  is  due  to  non-identical  laser  pulses  or  the  change  in  ambient 
conditions  as  a function  of  laser  interpulse  time.  Further  investigations 
would  be  required  to  determine  the  amplitude  sensitivity  on  laser  energy 
and  interpulse  time.  It  is  apparent  from  Eq.  (2.5)  that  two  key  parameters 
apart  from  irradiance  are  the  pressure  and  density  of  the  unshocked  gas. 

If  the  second  pulse  is  generated  before  the  perturbed  environm.ent  of  the 
first  pulse  has  returned  to  ambient  conditions,  a different  induced  surface 


pressure  time  history  will  result.  From  blast  wave  theory,  it  would  be 
expected  that  the  induced  surface  pressure  would  be  smaller  for  the  second 


pulse.  This  is  because  Q is  much  smaller  than  o and  p is  greater  than 


Pq.  However,  because  there  can  be  significant  mass  addition,  this  effect 


may  dominate  and  indeed  the  density  may  be  greater  than  o^.  Then  the 


second  pulse  really  is  traveling  into  a denser  gas  than  ambient  and  the 
induced  surface  pressure  could  well  be  enhanced.  This  behavior  will 
be  discussed  again  in  this  section  when  two  pulses  were  very  closely 
spaced. 

The  laser  system  assembled  to  investigate  multiple  pulse  radiation 
induced  sound  is  comp>osed  of  four  K-103  CO2  laser  kits  purchased  from 
Lumonics  Research  Ltd.  of  Canada.  These  four  TEA  lasers  are  charged 
from  one  high-voltage  power  supply,  and  triggered  individually  from  a 
four-channel  variable-delay  trigger  generator.  Therefore  one,  two,  three, 
or  four  laser  pulses  were  obtained  with  various  interpulse  times.  A 
schematic  diagram  of  the  experiment  appears  in  Fig.  3-1. 

The  CO^  TEA  laser  is  a short-pulse  laser.  Typically  the  output 
pulse  of  such  a laser  has  a very  short  high-power  spike  followed  by  a lower 
-power  but  longer  "tail".  This  "tail”  can  contain  up  to  2/3  of  the 
total  laser  energy.  The  experiments  consisted  of  operating  the  lasers  with 
a flowing  gas  mixture  of  He;N2:C02  of  5:1. 6:1.0  SCFH  measured  with 
flow  meters  calibrated  for  air.  The  laser  capacitors  were  charged  to 
35  kV,  and  had  an  output  of  approximately  7 J.  This  output  energy  was  the 
result  of  a tradeoff  with  pulse-to- pulse  stability.  The  pulse  shape  was 
monitored  with  an  L-401  pyroelectric  detector  manufactured  by  Barnes 
Engineering,  Figure  3-2  is  a representative  trace  of  the  detector  output. 
We  observe  a spike  and  a tail  of  several  usee  duration.  The  laser  energy 
was  monitored  with  a 50  D pyroelectric  energy  meter  manufactured  by 
Lumonics.  The  amplitude  of  the  approximately  50  msec  wdde  pulse  is 
proportional  to  the  laser  pulse  energy. 

The  four  channel  trigger  generator  was  designed  and  built  by 
Physical  Sciences  Inc.  A ''bounceless''  switch  utilizing  a monolithic  timing 
circuit  (Signetics  NE555  timer)  fed  four  556  dual  timers  whose  output 
signals  were  amplitifed  to  yield  the  required  trigger  voltage  amplitude 
(approximately  15  volts).  The  laser  beams  were  directed  to  focusing 
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mirrors  by  four  gold  coated  flat  mirrors.  Four  2.  5 m focal  length 
copper  mirrors  focused  the  laser  beams  to  approKimately  the  same  spot 
on  the  water  in  a large  tank.  The  worst  error  in  location  was  approxi- 
mately one-half  the  beam  size.  The  rectangular  laser  beam  at  the  water 

surface  was  0.  8 x 1.6  cm,  measured  from  a burn  pattern.  The  laser 

2 

energy  with  this  area  corresponds  to  a fluence  of  6 J/cm  . The  tank  was 
1.6  m in  diameter,  1.8  m high  and  it  was  filled  with  water  to  a depth  of 
1 . 4 m. 

The  acoustic  signature  was  detected  with  an  LC-5-2  hydrophone 

manufactured  by  Celesco,  The  hydrophone  has  a nominal  bandwidth  of 

500  kHz,  although  calibration  above  100  kHz  is  not  available  from  the 

manufacturer.  The  calibration  curve  is  reproduced  in  Fig.  3.  3.  From 

the  calibration  information,  the  average  free  field  sensitivity  is  -126.9  dB 

6 2, 

re  1 volt/ubar  which  corresponds  to  0.  452  V/10  dynes/cm  . The  cali- 
bration at  100  kHz  will  be  used  (-128  dB)  because  100  kHz  is  representative 
of  the  frequency  peak  of  the  spectrum.  The  sensitivity  is  0.  40  V/10^ 
dynes  / cm^. 

A series  of  one,  two,  three,  and  four  laser  pulse  experiments  was 
performed.  A number  of  representative  experimental  traces  will  be  dis- 
cussed in  detail.  Many  other  experimental  traces  were  obtained,  but  will 
not  be  discussed  in  this  report.  The  data  were  captured  by  photographing  a 
Tektronix  475  oscilloscope  trace.  The  resultant  photograph  was  digitized 
from  a magnified  copy  of  the  original  photograph.  A minimal  distortion  of 
approximately  1%  is  introduced  in  this  process.  The  waveforms  that  were 
analyzed  are  tabulated  in  Table  II  along  with  the  number  of  pulses  and  the 
interpulse  time.  The  amplitudes  will  be  evident  in  the  computer-generated 
plots  of  these  waveforms. 

The  waveforms  and  spectra  of  these  waveforms  will  now  be  dis- 
cussed. A single  pulse  waveform  is  shown  in  Fig.  3.4,  where  the  vertical 


-24- 


*.V. 


f 


\ 


Fiq,  3,3  Hydrophone  Calibration 


TABLE  II 


ANALYZED  WAVEFORMS 

2 

Laser  Fluence  Per  Pulse  = 6 J/cm 

2 

Threshold  Fluence  ~ 3J/cm 


1. 

Single  Pulse 

2. 

Double  Pulse 

At  = 15  Usee 

3. 

Double  Pulse 

At  = 20  |isec 

4. 

Double  Pulse 

At  = 45  Usee 

5. 

Double  Pulse 

^t  = 3 |isec 

6. 

Triple  Pulse 

At  = 45,  64  Usee 

7. 

T riple  Pulse 

At  = 157,  152  Usee 

8. 

Four  Pulse 

At  = 38,  33,  30  ^sec 

9. 

Four  Pulses  made 

equally  spaced.  At  = 31  ^sec 

4 

10. 

Four  Pulse 

At  = 115,  17.  5,  50  lise^ 

> 

11. 

Four  Pulses  made  equally  spaced.  At  = 17.5  Jisec 

( 

COMPUTER  SYNTHESIZED  WAVEFORMS 
USING  SINGLE  PULSE  SIGNATURE 

12. 

Single  Pulse 

( 

13. 

Double  Pulse 

At  = 31  pisec 

•« 

14. 

Triple  Pulse 

At  = 31  Usee 

15. 

Four  Pulse 

At  = 31  ^isec 
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scale  has  been  normalized  to  1.5  x 10’’  dynes /cm^.  The  experimental 
data  were  obtained  with  the  hydrophone  at  a depth  of  1 m.  The  attenuated 
signatures  correspond  to  additional  propagation  distances  of  500,  1000, 

1500  m from  the  observation  point.  Spherical  spreading  is  not  included  in 
the  normalization  constants.  The  energy  spectral  density  of  this  waveform 
is  shown  in  Fig.  3.  5.  We  note  that  the  peak  in  the  spectrum  is  at  100  kHz. 
This  is  to  be  compared  with  the  theoretical  prediction  of  100  kHz  with  a 
3 usee  pulse,  as  shown  in  Fig.  2.6.  The  experimentally  observed  peak 
corresponds  to  a laser  pulse  length  of  approximately  3 psec,  which  was 
observed,  see  Fig.  3.2. 

We  now  compare  the  amplitude  of  the  signal  with  the  theoretical 

prediction.  These  will  be  evaluated  at  5psec  into  the  signal,  because  the 

compressive  spike  is  probably  not  accurately  resolved  by  the  hydrophone. 

4 2 

From  Fig.  2.2  we  find  that  at  5 ^sec  the  amplitude  is  8 x 10  dynes /cm 
2 2 

for  3 J/cm  (the  threshold  fluence  of  3 J/cm  was  subtracted  from  the 

2 2 
laser  fluence  of  6 J/cm  ).  This  is  to  be  compared  with  4.5  x 1 0'^  dynes /cm 

from  Fig.  3.4.  Note  that  both  experimental  theoretical  values  are  for 

r - 1 meter  in  this  comparison. 

The  square  of  the  digitized  waveform  has  been  summed  to  yield  an 

acoustic  energy  density.  Using  the  calibration  at  100  kHz  of  -128  dB  re  IV 

8 2 

per  ubar,  we  find  that  the  acoustic  energy  density  is  1.5  x lO"  J/cm^. 

Assuming  a cos  0 dependence  for  a radiation  pattern,  we  find  that  the 

_4 

total  acoustic  energy  is  3 x 10  J.  This  yields  a ratio  of  acoustic  energy 
to  laser  energy  of  4.  2 x 10 

A series  of  double  pulse  signatures  was  measured.  The  analyzed 
signatures  had  interpulse  times  ranging  from  3 to  45  )jsec.  We  will  first 
discuss  the  signatures  with  interpulse  times  ranging  from  15  to  45  usee. 
When  the  inverse  of  the  interpulse  time  is  less  than  the  frequency  at  the 
peak  of  the  spectrum,  a cosinusoidal  modulation  of  the  single  pulse 
spectrum  is  expected  (Eq.  (3.  1)).  This  modulation  frequency  is  equal  to 
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the  inverse  of  the  interpulse  time.  As  the  interpulse  time  is  lengthened, 
the  modulation  frequency  is  lowered.  Three  double  pulse  signatures  and 
their  respective  spectra  have  been  selected  to  show  this  trend.  These 
are  shown  in  Figs.  3.  6 to  3.11.  These  signatures  are  in  sequence  of 
increasing  number  of  peaks  in  the  spectra  as  the  modulation  frequency 
dec  reases . 

The  remaining  double  pulse  signature  to  be  considered  is  a special 
case.  The  interpulse  time  is  approximately  3 ^isec,  which  is  very  small. 
The  signal  is  shown  in  Fig.  3.  12,  and  the  energy  spectral  density  is  given 
in  Fig.  3.  13.  There  is  no  cosine  modulation  of  the  spectrum,  and  while 
not  an  ideal  example,  this  case  shows  that  wehn  the  interpulse  time  gets 
short  enough,  no  modulation  of  the  spectrum  occurs.  If  the  pulses  had 
been  identical,  the  spectrum  would  have  narrowed  significantly.  Ideally, 
the  spectrum  would  continue  to  narrow  in  proportion  to  the  inverse  of  the 
number  of  pulses.  The  utility  of  this  type  of  waveform  may  be  limited 
because  of  the  high  attenuation  of  the  signature  with  distance. 

Two  three-pulse  waveforms  were  analyzed.  The  first  is  shown  in 
Fig.  3.  14,  The  pulses  are  almost  but  not  quite  evenly  spaced,  i.  e,  , 

157  nsec  and  152  p.sec,  and  because  of  this,  the  ideal  modulation  of  a 
single-pulse  spectrum  illustrated  in  Fig.  2,8  is  not  expected.  Neverthe- 
less, we  do  see  this  type  of  modulation  in  portions  of  the  spectrum  in 
Fig.  3.  15.  We  do  not  expect  to  see  the  1 to  9 ratio  of  "bump"  to  peak 
found  in  Fig.  2.  8,  because  the  amplitudes  of  the  three  pulses  were  not 
identical.  The  second  three-pulse  signature  is  shown  in  Fig.  3.  16.  This 
is  an  illustration  of  a signature  with  unequal  pulse  spacing.  No  simple 
modulation  of  the  spectrum  is  expected  or  observed  as  shown  in  Fig.  3.  17. 
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The  first  four  pulse  signature  to  be  analyzed  is  shown  in  Fig,  3.  18. 
Unfortunately,  the  pulses  were  equal  neither  in  amplitude  nor  in  spacing. 

The  spectrum  is  given  in  Fig.  3,  19.  We  note  the  modulation,  but  also  the 
lack  of  any  coherent  modulation  of  the  spectrum.  The  pulses  were  arti- 
ficially made  equally  spaced  in  a computer  synthesized  signature.  This  is 
shown  in  Fig.  3.  <20.  The  spacing  was  made  equal  to  31  |jsec.  We  note 
the  more  or  less  coherent  modulation  of  the  spectrum  in  Fig.  3.2  1.  A 
modulation  such  as  shown  in  Fig.  2.8  is  not  expected,  because  the  amplitude 
of  the  pulses  are  not  equal  - causing  the  bumps  in  between  the  spikes  in  the 
spectrum  to  be  larger  than  for  an  ideal  case. 

The  second  of  the  four  pulse  signature  is  shown  in  Fig.  3.22.  The 
spacing  of  the  pulses  is  not  uniform,  so  a ragged  spectrum  is  expected. 

The  spectrum  is  shown  in  Fig.  3.23.  Indeed  the  spectrum  does  not  show 
any  repetitive  structure.  The  broad  band  character  is  maintained  by  the 
unevenly  spaced  pulses.  If  the  pulses  had  been  evenly  spaced,  the  spectrum 
would  be  much  more  characteristic  of  Fig,  2,8  for  four  cycles  of  a sine  wave. 
An  evenly  spaced  signature  was  synthesized  from  Fig.  3.22,  This  was  done 
by  making  the  pulse  separation  equal  to  17.  5 usee,  which  is  the  spacing 
between  the  second  and  third  pulse.  The  signature  is  shown  in  Fig.  3.24, 
and  the  spectrum  is  shown  in  Fig.  3,25,  We  note  the  rather  significant 
narrowing  of  the  main  spectral  components,  and  also  the  two  double  bumps 
between  the  peaks  in  the  spectrum.  These  bumps  are  very  characteristic 
of  a four  pulse  spectrum;  e.  g.  , see  Fig,  2,8,  The  lack  of  data  for  evenly 
spaced  pulses  is  not  due  to  any  fundamental  physical  limitation.  When  these 
data  were  taken,  the  lasers  did  not  trigger  with  a stable  interpulse  time. 

In  an  effort  to  illustrate  the  behavior  of  the  signature  and  spectrum 
with  pulse  number,  the  following  computations  were  made.  The  single  pulse 
signature  of  Fig.  3.4  was  used  as  a single  pulse  signature,  and  repeated  to 
form  a two  -.three  -,  and  four  - pulse  signature  (the  little  bump  in  the 
signature  at  40  usee  was  eliminated).  The  computed  results  are  given  in 
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Figs,  3,26  to  3,33,  and  the  characteristic  modulations  of  approximately 

sin^  (N-rf/f  )/sin^  (nf/f  ),  Also,  the  peaks  in  the  spectrum  are  increasing 
“ 2 ° 

approximately  as  N , as  expected.  In  addition,  the  ratio  of  the  peak  to 

2 

the  small  bumps  between  the  spikes  is  approximately  N , These  wave- 
forms and  spectra  are  analogous  to  those  for  Fig,  2,  8,  Scaling  can  be 
established  from  these  calculations. 

We  now  discuss  the  spectral  narrowing  as  a function  of  the 
number  of  pulses.  The  measure  of  the  narrowing  which  will  be  used  is 
the  ratio  of  the  modulation  frequency,  f^,  (the  inverse  of  the  interpulse 
time),  to  the  width  of  the  highest  peak  in  the  spectrum  measured  at  half 
amplitude.  The  ratio  of  the  peak  frequency  to  the  width  will  be  used  for  the 
single-pulse  cases  which  will  allow  comparisons  of  waveforms  made  at 
different  interpulse  times.  The  waveforms  described  in  Table  2 were  analyzed 
from  the  energy  spectral  densities.  The  results  are  shown  in  Fig,  3,  34, 

The  o's  are  the  experimental  data.  The  triangles  are  the  two  four  pulse 
cases  that  were  synthesized  by  making  the  pulses  equally  spaced.  The  q's 
are  the  synthesized  waveforms  from  the  single-pulse  waveform.  We  see 
the  spectrum  narrowing  as  the  number  of  pulses  increases.  The  best 
three  pulse  data  point  was  28%  broader  than  an  ideal  case,  a significant 
demonstration  of  our  ability  to  control  the  spectrum.  The  four-pulse  data 
signatures  that  were  corrected  by  making  the  pulses  equally  spaced  exhibit 
even  better  results.  These  cases  are  3%  and  13%  broader  than  the 
synthesized  signal  with  equal  amplitude  pulses.  Therefore,  the  primary 
goal  of  the  study,  which  was  to  demonstrate  the  ability  to  control  the  energy 
spectral  density  with  the  number  of  laser  induced  acoustic  pulses,  has  been 
achieved. 

We  have  restricted  the  plotting  of  the  energy  spectral  density  to  the 
input  data.  In  Fig,  3,  35,  we  give  an  example  of  the  change  in  energy  spectra  1 
density  with  range  from  the  single-pulse  case  shown  in  Fig,  3,4,  We  note 
the  expected  trend  in  the  peak  of  the  spectrum  to  lower  frequencies  for 
greater  propagation  distances, 
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IV,  SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  WORK 

The  basic  experimental  goal  of  the  study,  the  ability  to  control 
the  energy  spectral  density  of  a laser-induced  acoustic  signature,  was 
accomplished.  Specific  illy,  it  has  been  demonstrated  that  as  the  number 
of  evenly  spaced  laser  pulses  increased,  the  width  of  the  fumlamental  or 
harmonic  peak  in  the  spectrum  is  decreased.  The  theoretical  performance 
was  not  achieved  because  the  laser-induced  pulses  were  not  identical  in 
amplitude,  shape, and  interpulse  time.  However,  an  experimentally  ohtaineri 
three-pulse  signature  was  analyzed  in  detail,  and  it  was  found  that  the 
spectrum  was  only  1.28  times  as  broad  as  the  theoretical  value. 

Because  there  is  only  a small  difference  between  the  theoretical  width  ind 
the  experimental  width,  the  basic  objective  of  the  experiments  has  been 
achieved. 

The  secondary  goal  of  the  experiments  was  to  obtain  signatures 
which  did  not  have  even  interpulse  times,  so  that  the  spectrum  would 
remain  broad.  It  is  not  possible  to  fully  demonstrate  this  concept  w ith  only 
‘‘our  pulses,  but  it  was  demonstrated  that  the  resulting  spectrum  was 
much  broader  than  one  with  four  evenly  spaced  pulses. 

The  experiments  discussed  in  this  report  have  demonstrated  the 
concept  of  spectral  narrowing  of  laser-induced  acoustic  signatures.  Quan- 
titative experiments  are  yet  to  be  done  to  verify  several  scientific  issues. 
First,  the  efficiency  of  converting  laser  energy  to  useful  acoustic  energy 
should  be  experimentally  determined.  A modest  theoretical  effort,  using 
the  existing  models  to  estimate  conversion  efficiencies  as  a function  of 
needed  acoustic  signatures,  should  be  carried  out.  The  relationship  between 
laser  parameters,  breakdown  level, and  the  resulting  acoustic  signature 
remains  to  be  established.  Finally,  the  optimum  IR  laser  wavelength  to 
maximize  the  conversion  efficiency  should  be  demonstrated. 
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Second,  comprehensive  experiments  need  to  be  done  to  quanti- 
tatively validate  the  relationship  between  the  various  laser  parameters 
(i.e.  number  of  pulses,  energy,  interpulse  time,  etc.)  and  the  resulting 
signature  parameters,  i.e.  amplitude  and  spectrum. 

Third,  the  acoustic  radiation  pattern  needs  to  be  studied  because 
many  applications  would  utilize  the  acoustic  radiation  at  shallow  angles. 
Because  of  the  broadband  nature  of  the  source,  the  signature  and  radia- 
tion pattern  will  be  dependent  on  source  size  and  laser  parameters.  Both 
the  ideal  angular  dependence,  on  a flat  surface,  and  the  nonideal  case  of 
a wavy  surface  require  experimental  and  theoretical  investigation. 

After  the  above  issues  have  been  addressed,  the  needs  of  numerous 
applications  can  be  matched  with  capabilities.  A limited  systems  study  has 
already  been  carried  out  to  establish  some  specific  needs. ^ 
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